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Abstract 

Vacuum polarization screening corrections to the ground state en- 
ergy of two-electron ions are calcualated in the range Z = 20 — 100. 
The calculations are carried out for a finite nucleus charge distribu- 
tion. 

1 Introduction 

It is well known that the dominant contribution to the vacuum polarization 
correction in a strong Coulomb field arises from the Uehling potential [0]. The 
remaining part of this correction is called by the Wichman-KroU contribution 
[0]. Calculations of the Uehling contribution cause no problem and have 
been done by many authors. Accurate calculations of the Wichman-KroU 
contribution were carried out in ^ 

Recent measurements of the Lamb shift in highly charged ions |^ have 
shown a necessity of calculation of vacuum polarization screening diagrams, 
i.e. diagrams with one vacuum loop and one electron-electron interaction 
(Fig.l). In the present paper an accurate calculation of these diagrams for 
the ground state of a two-electron ion is given. Results of the calculations 
are compared with a previous evaluation of this correction 

Relativistic units {h = c = 1) are used in the paper. 

2 Basic formulas 

The vacuum polarization screening diagrams are shown in Fig.l. The formal 
expressions for the energy level shift due to these diagrams can easily be 



1 



derived using the two-time Green function method [y]. The contribution of 
the diagrams shown in Fig la is 



AEa = J2i-lf J2 {{PaPb\IiO)\nb)^^{n\U^p\a) + 

P r ^o, 

{PaPb I / (0) I an)^ — {n \ U^p \ b) + 
{Pa I U^p I n)^ — {nPb | / (0) | ab) + 
{Pb I U^p I — {Pan | / (0) | ab)}. (1) 

Here Uyp is the vacuum polarization potentiah 

oo 



-oo 



A* 

/(a;, |x-y|) = ay^^^^ exp(i | u; 1 1 x - y |), (3) 

a'^ = {l,Oi) , a are the Dirac matrices; a, b are the Is states with a spin 

projection m = ±i; P is the permutation operator: G(uj, x, y) = ^''^•'fJ*''^! 

is the Coulomb Green function. 

The contribution of the diagrams shown in Fig lb is 

AE, = {-if {PaPb I U'yp I ab), (4) 



where 



oo 



?7yp(x,y) = J du J dzi J dz2- "^'^ 



-oo 



|x-zi| |y-Z2| 
xTriaf'Giuj, zi, Z2)a''G(cj, Z2, z^)) . (5) 

The contributions (1) and (4) are ultraviolet divergent. The most simple way 
to renormalize these contributions is to expand the vacuum loop in powers of 
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the external field (aZ). In this case the contribution of the diagrams with an 
odd number of vertices in the vacuum loop (with free-electron propagators) 
is equal to zero according to the Furry theorem. In this expansion, only 
the first nonzero term, called as the Uehling term, is infinite. The charge 
renormalization makes this term finite and its calculation causes no problem. 
The higher orders (in aZ) terms are finite. However, regularization is still 
needed in the second nonzero term due to the spurious gauge dependent piece 
of the light-by- light scattering contribution. As it was shown in [jlO|, [11], |], 
in the calculation of the vacuum polarization charge density, based on the 
partial wave expansion of the electron Green function, the spurious term does 
not contribute if the sum over the angular momentum quantum number k 
is restricted to a finite number of terms (|k| < K). We found that this rule 
is also correct in the case of the diagrams shown in Fig. lb (see Appendix 
A). Thus the Wichman-KroU contribution is calculated by summing up the 
partial differences between the full contribution and the Uehling term. 

3 Calculation 

3.1 Diagrams "a" 

Renormalized expression for the Uehling potential is well known 



Ui^niir) = -az£Jdr'Arrr'p{r')Jdt{l + —''^^ 



2 



1 



t2 



[exp {—2m\r — r'\t) — exp (— 2m(r + r')t)] 
Amrt ' 

where \e\Zp{r) is the density of the nucleus charge distribution ( / p{r)dr = 
1). It is interesting to note that with a good precision (~ 0.3% for Z = 80) 
the formula (6) can be replaced by a simpler formula (see Appendix B): 



oo 



KeMir) ^ V{r)- J dt{l + — )^^^exp(-2mrt) , (7) 
1 

where V{r) is the potential of an extended nucleus. The calculations were 
carried out for the Fermi model of the nuclear charge distribution using the 
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exact formula (6). To calculate the reduced Green function which appears in 



the formula (1) the B-spline method for the Dirac equation |12| was used. A 
change of the result, due to a one percent variation of the root-mean-square 
charge radius of the nucleus, was chosen as the uncertainty. 

The calculation of the Wichman-KroU contribution caused no problem 
too. Rotating the contour of the uj integration in the complex uj plane along 
the imaginary axis one get the following equation 

I, , oo oo oo 

2a f f f 1 

Uwxi^) = — H I I / / dyy"^ / dzz"^ rV{z) 

vr J J J max{x,y) 





2 

X 

i,k=l 



Y: Re{F:'{tuj,y,z){Gf{zu,y,z) - F:\tu,y,z))), (8) 

where Gjf are the radial components of the partial contributions to the 
bound-electron Green function and F*^ are the radial components of the par- 
tial contributions to the free-electron Green function. The homogeneously 
charged spherical shell model of the nucleus was used in the calculation of 
Uyyj^. In this case the Green function is expressed analytically in terms of the 
Whittaker, Bessel, and Hankel functions [|10|, |^. The reduced Green function 
which appears in (1) as well as the Is state wave function were calculated 
for the Fermi model. The k series in (8) converges rapidly with k, increasing, 
so we needed |k| < 5 to reach the relative accuracy not worse than 10~^. 

3.2 Diagrams "b" 

The calculation of AEi, was carried out in the same way as AE^. The corre- 
sponding expression for the Uehling operator is also well known: 

oo 



f/U(x, y) = cc — j dt{l + _)^^-^ exp(-2m | x - y | t) . (9) 

I X y I ^ ^ 

As before, the calculation of this term was made using the Fermi model of 
the nucleus charge distribution. The uncertainty of the results was estimated 
in the same way as for the Uehling part of the diagrams " a" . 

The calculation of the Wichman-Kroll contribution of the diagram "b" 
was done by summing the partial differences between the expression (5) and 
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the corresponding expression with the bound-electron Green functions re- 
placed by the free-electron Green functions. However, we encountered some 
obstacles trying to implement this scheme directly. The origin of these ob- 
stacles is that, in the direct numerical calculation, it is difficult to reach a full 
cancelation of some terms with large magnitude. To solve this problem we 
used the Furry theorem. Taking into account that according to this theorem 
only terms containing even powers of Z give nonzero contributions to the 
diagram shown in Fig lb, the Coulomb Green function was divided into two 
parts, each containing only odd or even powers of the nucleus charge. After 
that the expression Tr{GkGk) was replaced with Tr{Gf'^Gf'^ + Gf^''Gl''^''). 
Then the last trace was calculated analytically (see Appendix C). To control 
the accuracy of this procedure we calculated the correction due to interac- 
tion with an additional external field AV = — - to the first order vacuum 
polarization contribution for a point nucleus (Fig. 2). In this case the sum 
of all the corrections must be equal to the value dEyp/dZ where Eyp is 
the first order vacuum polarization contribution. The results of the test are 
given in the table 1. The calculation of the Wichman-Kroll contribution of 
the diagram shown in Fig. lb was carried out for a point nucleus. However, 
because of the smallness of this contribution (~ 0.002 eV for Z=90), the 
finite nuclear size correction can be neglected. 

The complete results of the calculation are presented in the table 2. The 
values of the root-mean-square charge radii were taken from |T^. For com- 



parison, the results of a previous calculation done in [§[ are listed in the last 
column of the table. 
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A 



The spurious gauge dependent piece of the hght-by-hght scattering contribu- 
tion for the diagram shown in Fig. lb is given by 



AE 



M 



(10) 



M 



du J f/xirfx2A(f)(xi)A('')(x2) J dy,dy2 



2tt Jcf 

xTr[a^FM{uj, xi, yi)V{yi)FM{uj, yi, xa) 

xaTM^uj, X2, y2)V(j2)FM{uj, y2, xi) 

+a''FMiuj, xi, X2)a''FM{uJ, X2, yi) 

xV{yi)FM{uj, yi, y2)V^(y2)i^Af(c^, y2, xi) 

+a''FM{uj, xi, yi)V{yi) FMiuj, yi, y2) 

x\/(y2)FAf(u;, y2, X2)a''FM(u;, X2, xi)]. (11) 

J dy^l^{y)a'^D,^{0,^-y)My), D,^ {u, x-y) is the photon 
propagator, Fm{uj, x, y) is the free electron propagator with the electron mass 
replaced by a hypothetical heavy mass M. Calculating a limit as M becomes 
infinite yields 



where ylff ^ fx 



lim 

A/^oo 



^-J^ du J d^,Ai^\^i)A^^\^,)V\^^] 



Ani Jcf 
xTr 



lim 

X2^Xl 



du^ 



F(u;,Xi,X2) 



e = ^j^^dujd^^J^;\^MVi?^i)V\^i)\}rnJdy 



xTr 



a 



duj'- 



:F(a;,Xi,y)a'=F(^,y,X2^ 



(12) 



(13) 



where F(c(j,x, y) is the free-electron propagator; i,k = 1,2,3. The term ^"-^ 
can be treated in the same way as the light-by-light scattering graph in the 
third-order vacuum polarization [|l^, So, we restrict our consideration to 
the term C,^. Using an identity = i[HF,x''], where Hp is the free-electron 
Hamiltonian, we find 



— / du 

Arc Jcf 



J d^,At\^i)A'i\^,)V\^i) hmy. 
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xrr(a^— F(a;,Xi,X2)) 



= i- / dxiAf Hx04'^(xi)y^(xO lim lim (xt - 4) 

47r J 12— >oo X2— »xi 



X rr(a'^F(a;,xi,X2)) 



(14) 



uj=—iCl 

The first derivative of F with respect to u; is 

r(i ^ (i iexp(— , , 

-ia-x + Pm + uj-\ (15) 

LUX-' Att a 



duj 



where x = Xi — X2, x = \x.\, d = \Jvn? — cu^. Wc can sec that is equal 
to zero independently of the order of taking limits in (A5). In the partial 
expansion the terms F'^iuj, x, y) are finite together with their first derivatives 
with respect to uj. It follows that the sum of a finite number of terms in the 
K series is equal zero. 

B 

Let us represent the Uehling potential in the form: 



UueHiir) = -aZ-ydt— ^(1 + — )^^^exp(-2mrt)A(r,t),(16) 

1 

where 

/•CO 

A{r,t) = / dr'r' p{r')[exp {—2m{\r — r'\ — r)t) — exp (— 2mr't)]. (17) 
Jo 

Expansion of the exponents in the Taylor series yields: 

[exp {-2m{\r - r'\ - r)t) - exp {-2mr't)\ = Umr< + 0{m^r'r<f). (18) 
Taking into account that 

/ dr'p{r')r'r< = -Vir) (19) 

Jo 4:7raZ 

one can get 

U^eHiir) ^ V{r)- j dt{l + — )^^^exp(-2mrt) . (20) 
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c 



The analytical expressions for the radial components of the Coulomb Green 
function are well known (see , e.g., [14]). In the units H — c — m — 1 for 
xi < X2 we have {G = {lu — H)~^) 



Gl}{u, xi,X2) = -{l + u)Q 



(A - iy)M,_.^,{2dx,) -{k- —)M,^.^,{2dx, 



X 



aZ 



+ -j^)W,_i^^{2dx2) + W^,+ i,,(2dx2) 



{u,Xi,X2) 



-dQ 



aZ , 



(A - iy)M,_i j2dx,) -{k- —)M,^i A2dxi) 



X 



aZ , 



[K+—)W,_^^^{2dx2) - W,^.^,{2dx2) 



G^{u;, xi,X2) = -dQ 



X 



ckZ 

(A - u)M,_i^,{2dx,) + {k- — )M,+ i^,(2cixi) 
aZ , 



[K+—)W,_.^^{2dX2) + W,^^^{2dX2) 



Gi'{u;,x^,X2) = -{l-u;)Q 



otZ 

(A - v)M,_._^^{2dx,) + {k- —)M,^^,{2dxi) 



X 



aZ 



K+—)W,_^^^{2dX2) - W,^.^^{2dX2) 



(21) 



where d = a/I — 0;^, A = {k^ — (aZ)'^), v 



aZui 



, Q 



Ma,pi and Wa,p are the Whittaker functions. For xi > X2 the radial Green 
functions can be obtained from the symmetry condition 



Defining 



Gf{uJ,Xi,X2) = G''^\uJ,X2,Xi). 

A = -Q{X-iy)M,_r^,{2dxr)W,_^^,{2dx2), 

B = -Q{\-v)M^_._^^{2dx^)W,^^y^{2dx^), 

C = -QM,+i ,(2rf,xi)l^,_i,,(2rfa;2) , 

D = -QM^^U2dx,)W,^.A2dx2), 



(22) 



and taking into account that for uj = ie, where e is real, changing the sign 
of Z is equal to replacing A, B, C, and D with their complex conjugated 
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values, we find for the radial parts of the Coulomb Green function containing 
only odd or even powers of Z the following expressions 



G 



11, odd 



G 



12,odd 



G 



21,odd 



22,odd 



/^12,even 



/^2l,even 
^ K 

Q22,even 



1 + ^){%{Klm{A -D) + Im{B) - -fIm{C)) 
"^^ReiA + D)}, 



d 



aZ 



d{iUlm(A + D)- ImiB) - ^ImiC)) + —ReiA - D)\ , 

a 

olZ 

d{i{Klm{A + D) + Im{B) + -fIm{C)) + —rRe{A - D)] , 

(Jj 

(1 - Lu){i{Klm{A -D)- Im{B) + 'fIm{C)) 
+^Re{A + D)}, 



d 



(23) 



^ii,e^en ^ {1 + u;){KRe{A - D) + Re{B) - jRe{C) 



+i^Im{A + D)}, 
d 



aZ 



d{KRe{A + D)- Re{B) - ^Re{C) + i—Im{A - D)} , 

(Jj 

OlZ 

d{KRe{A + D) + Re{B) + -fRe{C) + i—Im{A - D)} , 



= {1 - u){KRe{A - D) - Re{B) + -fRe{C) 

.aZ ^ . , 
+t—Im{A + D)}, 



(24) 



where 7 = — ^"J^ . 

In the calculation of the Wichman-Kroll part of the diagram shown in 
Fig. lb the expression Re{Tr{a^G{uj ,x,y)a^G{uj ,y,x))) appears. After hav- 
ing integrated over the angles we have to evaluate the following expressions 



^0 = Re{ E E(G^^'('^,^,1/))'), (25) 

sign{K)=—l i,fe=l 

^1 = Re{ E (G;.^G- + G-G^^ + G-G- + G-G-)) (26) 

sign{K)=—\ 

S2 = 2Re( E (Gl'G'J + GfGD + Gl'GlUGfG'J)), (27) 

sign{K)=—l 



9 



where k' = —sign{K,){\K,\ + 1). Using (C3)-(C4) and the Furry theorem we 
transformed these expressions to the following ones 

50 = 8{{^^^ + K''){Re{Af + Re(Df - Im{Af - Im{Df) 

+Re{Bf - Im{Bf + -i{Re{Cf - Im{Cf)} 

-16£^^{i?e(A + D){Im{B) - -fIm{C)) 
+Im{A + D){Re{B) - 7i?e(C))} 

-16£'{(^^|£ + K^){Re{A)Re{D) + Im{A)Im{D)) 
+-f{Re{B)Re{C) + Im{B)Im{C))} (28) 

51 = 8d''{{K''+^^^){Re{Af + Re{Df-Im{Af-Im{Df) 

+7(7m(C)' - ReiC)^) + Im{Bf - Re{Bf} (29) 

S, = IM\{\k\{\k\ + 1)-^^) 

x{Im{A)Im{A') + Im{D)Im{D') - Re{A)Re{A') - Re{D)Re{D')) 
+Im{B)Im{B') - Re{B)Re{B') + 77' (/m(C)/m(C") 
-Re{C)Re{C'))} (30) 

The equations (28)- (30) were used in the numerical calculation. 
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Table 1: Control calculation of the correction due to interaction with an ad- 
ditional external field AV{r) = —a/r to the first order vacuum polarization 
contribution for the Is state (see Fig. 2) (in eV). The label "a" corresponds 
to the sum of the first two diagrams in Fig 2. The label "b" corresponds to 
the third diagram in Fig 2. The calculation is carried out for a point nucleus . 



z 










AEvP 


cH^Jyp 
dZ 


20 


-0.019183 


-0.006487 


0.000072 


0.000076 


0.025522 


0.025521 


40 


-0.16181 


-0.05215 


0.00206 


0.00221 


0.20970 


0.20969 


60 


-0.67136 


-0.19446 


0.01650 


0.01702 


0.83231 


0.83234 


90 


-4.5280 


-0.9586 


0.1970 


0.1711 


5.1185 


5.1186 


100 


-9.2613 


-1.6460 


0.4532 


0.3543 


10.1000 


10.0999 



Table 2: Contribution of the vacuum polarization screening diagrams to the 

ground state energy of two-electron ions (in eV). 



z 






^El,, 


^E^K 


^E^K 


AEyp 


Ref. [8] 


20 


3.478 


0.0090 


0.0010 


-0.0000 


-0.0000 


0.0100 




30 


3.928 


0.0316 


0.0035 


-0.0003 


-0.0000 


0.0348 




32 


4.072 


0.0388 


0.0043 


-0.0004 


-0.0000 


0.0427 


0.0 


40 


4.270 


0.0807 


0.0091 


-0.0011 


-0.0000 


0.0887 




50 


4.655 


0.176 


0.020 


-0.003 


-0.000 


0.192 




54 


4.787 


0.234 


0.027 


-0.005 


-0.000 


0.255 


0.2 


60 


4.914 


0.350 


0.040 


-0.009 


-0.000 


0.380 




66 


5.224 


0.515 


0.058 


-0.016 


-0.000 


0.557 


0.6 


70 


5.317 


0.661 


0.074 


-0.022 


-0.000 


0.713 




74 


5.373 


0.845 


0.093 


-0.030 


-0.000 


0.908 


0.9 


80 


5.467 


1.215 


0.132 


-0.049 


0.000 


1.298 




83 


5.533 


1.455 


0.156 


-0.062 


0.001 


1.550 


1.6 


90 


5.645 


2.214(1) 


0.230 


-0.105 


0.002 


2.341(1) 




92 


5.860 


2.493(2) 


0.256 


-0.122 


0.003 


2.630(2) 


2.6 


100 


5.886 


1.063(1) 


0.398 


-0.221 


0.008 


1.218(1) 





Figure 1: Vacuum polarization screening diagrams 













o 



<3 



J\A/VX 



JXA/VX 



b 



Figure 2: Diagram equation for control calculation of the correction due to 
interaction with an external field AV — —a/r to the first order vacuum 
polarization contribution. 




